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分 段 复合 多 尺度 模糊 炉 和 IGWO-SVM 的 脑 电 情感 识别 


C 吴 清 
(河北 工业 大 学 人 工 智 能 与 数据 科学 学 院 ， 天 津 300401) 


摘 要 : 脑 电 情感 识别 在 辅助 测 说 、 临 床 医 学 、 人 机 交互 中 发 挥 着 重要 作用 。 为 提高 脑 电 的 情感 识别 率 ， 提 出 了 分 段 
复合 多 尺度 模糊 蛟 算法， 采用 分 段 粗 粒 化 和 计算 复合 多 尺度 模糊 蚁 的 策略 ， 使 提取 的 特征 较 好 的 解决 了 数据 缺失 和 计 
算 不 准确 的 问题 ; 同时 构造 了 应 用 余弦 非 线 性 收敛 因子 和 动静 态 位 置 更 新 的 灰 狼 算法 优化 支持 向 量 机 的 分 类 模型 。 为 
证 明 所 提 的 两 种 算法 的 有 效 性 ， 进 行 了 仿真 实验 的 验证 ， 并 在 公开 DEAP 数据 库 下 与 几 种 常见 的 支持 向 量 机 优化 模型 
比较 脑 电 的 情感 识别 率 , 结果 表明 在 提出 的 模型 下 , 效 价 、 唤 醒 度 、 优 势 度 、 喜 欢度 的 平均 识别 率 分 别 为 87.27%、87.81%、 
89.06%、87.58%， 均 高 于 其 他 算法 。 另 外 对 比 了 高 / 低 喜 欢度 下 效 价 和 唤醒 度 的 分 类 ， 实验 说 明 喜 欢度 低 时 情感 识别 率 


较 高 。 
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EEG emotion recognition based on piecewise complex multi-scale 
fuzzy entropy and IGWO-SVM algorithm 


Wei Xue, Wu Qing 
(School of Artificial Intelligence & Data Science Hebei University of Technology, Tianjin 300401, China) 


Abstract: EEG(Electroencephalography) emotion recognition are important in auxiliary lie detection, clinical medicine and 


human-computer interaction. The paper used a piecewise complex multi-scale fuzzy entropy to increase the recognition rate. It 


pl. proposed piecewise coarse graining and calculating complex multi-scale fuzzy entropy. It could solve the problem of missing 


data and calculating inaccurate entropy. This paper constructed an IGWO-SVM classification model which used cosine 


nonlinear convergence factor and dynamic and static position. To prove the validity of the algorithms, The paper performed the 
simulation experiment. Comparing the EEG emotion recognition rate with several common SVM optimization models in the 
public DEAP database, the results show that the average recognition rates of valance, arousal, dominance and liking are 
87.27%, 87.81%, 89.06% and 87.58% in proposed model, which is higher than other algorithms. Comparing the classification 
of valance and arousal under high/low liking, experiments show that the recognition rate is higher when the liking is lower. 

Key words: EEG signals; emotion recognition; improved grey wolf optimizer algorithm; SVM optimization algorithm; 


piecewise complex multi-scale fuzzy entropy 


A, Tees REE AN l 2007 年 Chen EAE SPRANG 

的 改 ; 模糊 箭 算法 辐 ， 这 是 一 种 新 的 序列 复杂 度 的 测量 方 

情感 识别 是 通过 对 传感器 采集 的 信号 进行 分 析 处 理 从 而 得 法 ， 并 成 功 地 运用 在 了 体 表 肌 电信 号 的 特征 提取 和 分 类 。 本 文 

出 被 采集 者 的 情感 状态 的 过 程 。 基 于 生理 信号 的 情感 识别 方法 TE SAREREA fe h AL RERA A (multi-scale fuzzy 
比 其 他 指标 如 面部 表情 、 身 体 姿 势 、 语 音 语 调 等 更 难以 伪装 ， entropy, MSFE ) , KJE T DREE L RE BUNK BIE (piecewise 
也 能 提供 更 丰富 的 信息 。 基 于 脑 电 生理 信号 的 情感 识别 已 经 广 complex multi-scale fuzzy entropy, PCMSFE) 作为 脑 电 信和 号 的 
泛 用 于 辅助 测度 中、 心理 疾病 的 识别 所 、 与 由 金森 病人 间 的 人 IIE HEH T OBRAR S E RERA REINE. 
机 交互 BI 等 方面 。 支持 向 量 机 (support vector machine, SVM) 是 脑 电 信号 
近年 来 ， 基 于 非 线 性 动力 学 的 特征 提取 方法 在 脑 电 等 生理 情感 识别 中 应 用 最 多 的 一 种 方法 ， 它 可 以 解决 小 样本 、 非 线性 
言 号 的 研究 中 取得 了 较 好 成 果 ， 如 分 形 维 数 申 、 排 列 粹 申 、 谱 问题 ， 并 且 在 高 维 的 模式 识别 中 起 到 了 重要 作用 ,但 SVM 的 
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内 部 参数 对 其 分 类 效果 影响 很 大 [由 。 因 此 ， 很 多 研究 者 将 布谷 
鸟 搜索 算法 (cuckoo search, CS) 021、 粕 子 群 优化 算法 (particle 
swarm optimization, PSO )03、 遗 传 算法 (genetic algorithm , GA) 
U4] 324) HELL (differential evolution, DE) 051 等 算法 应 用 到 SVM 
的 参数 寻 优 上 ， 可 是 这 些 算法 对 于 SVM 参数 的 优化 效果 仍 有 
其 适应 性 的 问题 。Mirjalili 等 人 051 创 建 的 灰 狼 优化 算法 (grey 
wolf optimizer, GWO) 在 寻 优 方 面具 有 原理 简单 、 需 要 调整 的 
参数 少 、 易 于 实现 等 优点 ， 然 而 GWO 算法 也 存在 过 早 收敛 、 
容易 陷入 局 部 最 优 、 对 初始 种 群 有 依赖 性 等 问题 。 为 解决 上 述 
问题 ，Zhu 等 人 0 选取 了 一 种 DE 和 GWO 相 结合 的 混合 算法 ， 
不 过 两 种 算法 的 融合 会 导致 算法 体系 过 于 庞大 且 寻 优 效果 不 大 
里 想 。Saremi 08 将 动态 进化 种 群 策 略 应 用 到 GWO 算法 中 ， 提 
高 GWO 算法 的 全 局 搜索 能 力 。 因 此 ， 本 文 创 建 了 一 种 基于 余 
弦 非 线性 收敛 因子 和 动静 态 结合 进行 位 置 更 新 的 改进 灰 狼 算法 
(improved GWO, IGWO) ， 来 平衡 算法 的 局 部 搜索 能 力 和 全 
局 搜索 能 
本 文 主要 针对 脑 电 信号 情感 识别 中 涉及 的 特征 提取 和 分 类 
问题 进行 了 如 下 工作 : a) 为 防止 数据 的 缺失 和 计算 的 粒 不 准确 
的 问题 的 出 现 ， 采 取 了 一 种 分 段 粗 粒 化 和 计算 复合 粒 的 多 尺度 
WAM SE, Xt MSFE #0 PCMSFE 进行 了 白 噪 声 和 1 噪声 的 仿 
真 测试 和 脑 电 信号 的 情感 识别 测试 , 实验 证 明 PCMSFE 的 识别 
准确 率 高 于 MSFE;b) 对 于 脑 电 特征 的 分 类 ， 本 文 构建 了 一 种 
改进 的 灰 狼 算法 优化 SVM 的 参数 模型 ， 利 用 余弦 非 线性 收敛 


因子 和 动静 态 结 合 的 位 置 更 新 策略 进行 改良 ,选择 9 个 基准 测 
试 函数 进行 了 仿真 实验 ， 并 且 对 比 了 在 PCMSFE F, IGWO 和 


GWO, PSO, FA (firefly algorithm， 莹 火 虫 算法 ) 、GA 对 SVM 
优化 的 情感 识别 准确 率 ， 结 果 表 明 IGWO 对 SVM 进行 优化 的 


编号 ， 以 通知 被 试 者 实验 进 
乐 视 频 的 显示 ;d) 被 试 者 进 


实验 包括 以 下 几 个 步骤 ， 
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AN 


FRAT 


。 实 验 前 每 个 被 试 者 都 会 被 问 及 是 否 舒 适 ， 并 在 必要 时 进行 


a) 在 屏幕 上 显示 2s 的 当前 视频 


度 ;b) $s 的 基 
行 唤醒 度 、 效 价 、 


线 记 录 :;c) 1 min 的 音 


喜欢 度 和 优势 度 


的 自我 评估 。 经 过 20 次 实验 , 被 试 者 进行 休息 , 实验 者 对 信号 


的 质量 和 电极 的 位 置 进行 检查 。 
实验 数据 集 已 经 经 过 了 简单 的 预 处 理 ， 将 数据 下 采样 到 
128 Hz， 并 进行 了 4.0~45 Hz 的 滤波 处 理 ， 使 用 盲 源 分 离 技 术 
去 除 眼 电 伪 迹 。 每 组 数据 记录 总 共 63 s， 包 括 60 s 实验 过 程 数 
据 和 3 s 的 基线 数据 。 本 文采 用 60 s 的 数据 段 进行 实验 。 
1.2 脑 电信 号 处 理 过 程 
般 的 脑 据 的 预 处 理 、 数 


括 了 对 数据 的 分 割 和 僵 加 ， 


信号 情感 识别 过 程 主要 包括 数 
类 


据 特征 的 提取 、 数 据 特征 的 分 类 


本 文 对 脑 电 


H 号 的 处 里 还 包 


以 提高 特征 数量 ， 


避免 因数 据 量 过 


小 而 产生 分 类 效果 不 佳 等 问题 。 本 文 对 脑 电信 号 的 主要 处 理 步 


又 和 具体 流程 框架 如 图 1 所 示 。 
V 
图 1 本 文 脑 电 信号 情感 识别 的 流程 框架 
a) 数据 预 处 理 。 脑 电信 号 的 噪声 源头 可 能 是 由 周围 设备 
产生 的 静电 或 电磁 场 造成 的 ， 除 了 这 些 外 部 噪声 外 ， 脑 电信 和 号 
还 受 身 体 运动 和 用 眼 动 作 等 人 为 因素 的 严重 影响 中 。 对 此 ， 
DEAP 数据 库 已 经 进行 了 一 些 相应 的 处 理 。 根 据 以 往 的 经 验 ， 


识别 准确 率 高 于 其 他 算法 。 另 外 ， 比 较 了 分 别 在 喜欢 度 高 和 喜 
欢度 低 的 情况 下 ， 效 价 和 唤醒 度 的 分 类 情况 ， 实 验证 明 喜 欢度 
低 时 效 价 和 唤醒 度 的 识别 率 更 高 。 

1 ”相关 算法 


1.1 


DEAP 数据 集 
DEAP (database for emotion analysis using physiological 
signals) 数据 库 09 是 一 个 公开 数据 库 ， 是 由 Koelstra 等 人 提出 
的 对 人 类 情绪 状态 分 析 的 多 通道 数据 集 。DEAP 数据 库 集 脑 电 
言 号 、 生 理 信 号 和 视频 信号 于 一 体 ， 提 供 公开 的 情感 分 析 数 据 
Se, 供 人 们 进行 实验 。 数 据 库 记录 了 32 名 被 试 者 的 脑 电 图 和 外 
恒生 理 信号 。 
32 位 健康 被 试 者 参与 了 实验 , 被 试 者 在 实验 之 前 被 告知 实 
验方 案 和 自我 评估 表 的 含义 。 被 试 者 在 实验 中 观看 40 个 能 够 引 
起 情绪 显著 变化 的 音乐 视频 。 实 验 在 有 可 控 照 明 的 两 个 实验 室 
环境 中 进行 ， 使 用 Biosemi ActiveTwo 系统 在 专用 记录 PC 
(Pentium 4, 3.2 GHz) 上 记录 脑 电信 号 和 外 围 生理 信号 。 音 乐 
视频 在 17 英寸 的 屏幕 上 呈现 ， 为 了 最 大 限度 的 减少 眼球 运动 ， 
所 有 视频 刺激 大 致 填 满 2/3 的 屏幕 ， 被 试 者 大 约 在 离 屏 幕 1 米 


在 实验 中 Alpha、Beta 波 对 结果 的 影响 较 大 ， 


较 小 20， 所 以 只 提取 这 两 个 


波段 进行 实验 。 通 


其 他 波段 的 影响 
过 带 通 FIR 滤波 


器 对 信号 进行 过 滤 ， 去 除 小 于 8 Hz 和 大 于 30 Hz 的 噪声 。 
b) 数据 分 割 。 分 别 对 Alpha, Beta 波段 的 数据 进行 分 割 ， 


选取 20s 的 移动 窗口 对 60 s 
共有 五 段 数据 ， 每 段 数据 为 
c) 特征 提取 。 选 择 脑 


(4 对) FP1-FP2, AF3-AF4, 


段 数据 中 8 个 通道 的 功率 谱 


数据 段 进 行 分 割 ， 
20 S。 


一 | 


HH 10 s, 


BIA J x 


RPE. 并 计算 MSFE 


特征 ， 以 供 后 续 比 较 不 同 特 


d) 数据 县 加 。 把 每 段 数 据 提取 的 12 MEME Te, JE 


征 对 识别 率 的 影响 


每 段 数据 并 行 合 加 ， 每 组 数 


个 Beta 波段 特征 =24 “EFI 


*5=200 个 样本 ,随机 
进行 测试 。 


据 


显著 影响 的 8 个 通道 
F3-F4, F7-F8 提取 特征 ， 提 取 每 


特征 和 PCMSFE 


o 


得 到 12 个 Alpha 波段 特征 +12 
E， 每 个 被 试 者 有 40 (40 个 实验 ) 


由 取 其 中 80% 的 样本 进行 训练 ,剩余 的 20% 


e) 特征 分 类 。 运 用 IGWO 优化 的 SVM 参数 模型 (定义 


为 IGWO-SVM 模型 ) 对 脑 


信号 进行 情感 识 


别 ， 对 数据 库 的 


效 价 、 唤 醒 度 、 优 势 度 、 喜 欢度 进行 分 类 , 将 GWO、PSO、FA、 
GA 优化 的 SVM 模型 与 IGWO-SVM 模型 进行 


对 比 。 同 时 ， 比 
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较 IGWO-SVM 模型 分 别 在 MSFE 特征 和 PCMSFE 特征 下 的 情 


ax iver ERB TH 


w 雪 ， 等 : 分 段 复合 sranmg RXR te 感 识别 


RA 


复合 多 尺度 模糊 焙 的 计算 办 法 。 分 段 粗 粒 化 可 以 在 数据 长 度 不 


感 识别 率 。 另 外 ， 利 用 IGWO-SVM 模型 比较 喜欢 度 高 低 不 同 
时 效 价 和 唤醒 度 的 分 类 效果 。 
1.3 AREASRERMM 

电信 号 是 非 稳定 的 随机 信号 ， 使 用 传统 的 波形 、 


地 
六 


是 尺度 因子 整数 倍 时 ， 解 决 数据 流失 问题 ; Le RRR 
算 能 够 解决 短 时 间 序 列 在 尺度 因子 大 时 计算 不 准确 的 问题 ， 并 
且 降 低 诱发 未 定义 灶 的 概率 。 主 要 处 理 过 程 如 下 

a) 原 始 序列 分 段 粗 粒 化 。 传 统 粗 粒 化 将 原始 时 间 序 列 按 尺 


7 


析 很 难 达到 人 们 预计 的 要 求 。 因 而 ， 非 线性 动力 学 特征 
到 脑 电信 号 的 研究 中 。 灶 是 度量 大 脑 复 杂 度 的 重 
WE AVL IRR ALVA ASB ae EY Ti) BSR PETIA, RPP AE 
Bea. HBOR. ELT DAI MT BSB BIE SB — Be ; 
BEAST ELIN iS — TP ADL, 否则 会 使 计算 没有 意义 
TRI RE RERA E, ASAT DA EBS I EE 
1, GAA BRACE. ATRAE ARE EE 
ENEEIER ERTA, KELOR h TERERAA 
的 概念 ， 用 来 衡量 不 同 尺 度 因 子 下 时 间 序 列 的 复杂 性 。 可 这 种 
算法 容易 导致 数据 丢失 并 使 计算 的 箭 不 够 准确 ， 因 此 本 文选 用 
了 一 种 基于 分 段 复 合 策略 的 多 尺度 模糊 烂 算法 进行 脑 电 信号 的 
特征 提取 。 

1.3.1 AM SR LARA 


te RB 
Ho 


WE 
ys 
洋 


m 


MSFE 方法 首先 将 原始 序列 {% :1<i<N} 粗 粒 化 ， 建 立新 


的 粗 粒 向 量 y,(z) 。 


T i=( 六 JrtHl 
其 中 : 5 =1,2,.…. 是 尺度 因子 。 当 7=1 时 ，y,(1) 是 原 序 列 。 
a rs 
将 粗 粒 序列 顺序 排列 后 ， 减 去 序列 均值 ， 形 成 向 量 Y”。 


PARERA (dr nr) 定义 矢量 Y” 和 YY 的 相似 度 Dy。 给 定 


PR RCA KA BER m 和 相似 容 限 r， 定 义 函 数 


m 1 N-m 1 N-m 
d ae Nomi 
jai 


” (2) 


类 似 地 ， 可 扩展 到 维 数 m+, 


1 


办 (m7)= N- “Sits m— Som (3) 


EME REO 


MSFE(m, n, r) = lim [In(9" (n, r)) 一 In (g (n, r))] (4) 


当 N 为 有 限 数 时 ， 式 〈4) 可 表示 为 
MSFE(m,n,r,N) = In(g” (n, r)) 一 In( 乡 ” (n, r)) (5) 
MSFE WAR RERIP PHA, CRT AIRE 


子 下 时 间 序 列 的 复杂 性 和 自 相 似 性 。 但 式 “5)〉 定义 的 MSFE 
容易 导致 数据 缺失 、 计算 结果 不 准确 , 或 导致 计算 的 炉 未 定义 。 


1.3.2 基于 分 段 复合 策略 的 多 尺度 模糊 炳 
为 解决 上 述 MSFE 存在 的 问题 ， 本 文采 用 了 分 段 粗 粒 化 和 


度 因 子 压缩 时 ， 如 果 数 据 长 度 不 是 尺度 因子 的 整数 倍 ， 就 会 导 
致 部 分 数据 的 缺失 。 受 多 斥 度 样本 精 的 粗 粒 化 改进 的 启发 己 ， 
将 分 段 粗 粒 化 方法 用 于 多 尺度 模糊 焙 ， 定 义 


ae 1 ie 
y? ie) Xp. 
T k-(j-1)r++p 


1s7<| E|, p= 01,2001 
T 


(6) 


其 中 :7 是 尺度 因子 ,WN 为 时 间 序 列 的 长 度 。 将 粗 粒 化 后 的 序列 


y(i) 按 顺序 排列 ， 计 算 向 量 Y”， 


y(i):1<i<|(N - p)/t |, p =0,1,2,..., -1 (7) 
=) mm) (8) 
i=1,2,..,.N—m+l1 
Y” 代表 从 第 i 个 点 开始 的 连续 m 个 y 的 值 减 去 均值 Yoli) 
)=2 Soi) 
AY” 5Y" 间 的 距离 4| 7”,7” | 为 两 个 对 应 元 素 的 最 大 差 
值 ， 即 
dy =d[¥".¥" |= max, {| y(i+k)—yo@O—y(F+k)— yo DI} 而 


i,j =1,2,..N—m,i¢ j 


H 


使 用 模糊 函数 (dn) TREY,” 和 了 ”的 相似 度 DY 


pr =p(dr mr)=e ts) 


y 


(10) 


其 中 :2 和 了 分 别 为 指数 函数 fa dy yn, r) 边界 的 梯度 和 宽度 。 


b)iTHA ZS RRO. MSFE 在 尺度 因子 过 大 时 会 导 
致 数据 长 度 过 小 ， 因 此 ，MSEFE 分 析 短 时 间 序 列 时 可 能 存在 计 
算 结 果 不 准 确 的 问题 。 虽 然 利 用 文献 [23] 中 的 方法 可 使 问题 得 
到 改善 ， 但 仍然 存在 计算 不 准确 和 诱发 未 定义 烂 的 可 能 性 。 本 
Mish THRE L RE BOM TIS CMSFE (complex 
multi-scale fuzzy entropy) 来 提高 计算 精度 ， 其 定义 为 


pak) 


CMSFE( y,z,m, r) = Ing (n, r)) -In(g" (n, r)) (11) 
h (n.r) r =e" (n.r) r 
oy" (nr)= BX + (n,r) 


p: g" (nr) AO" (mr) 可 分 别 通过 式 (2) (3) 计算 得 到 。 


式 (11) 可 进一步 被 简化 为 


:201808.00088v1 
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CMSFE(y,t,m,r)= "(P(e (n.r))- inf (6 “Oo (12) 


式 (12) 是 将 传统 MSFE 相 加 平均 得 到 的 ， 这 样 可 以 降低 
产生 不 准确 烂 和 未 定义 焙 的 概率 。 
1.4 灰 狼 算法 

SVM 在 脑 电 信号 的 特征 分 类 中 应 用 广泛 ， 由 于 SVM 分 类 
效果 受 内 部 参数 影响 很 大 。 相 对 于 PSO, GA 算法 , 使 用 GWO 
优化 SVM 需要 调整 的 参数 少 、 原 理 简 单 。 但 GWO 也 存在 容 
易 陷入 局 部 最 优 、 位 置 更 新 不 灵活 的 问题 ， 于 是 本 文 利用 余弦 
非 线 性 收敛 因子 来 平衡 全 局 搜索 和 局 部 搜索 ， 选 取 调 和 平均 数 
和 静态 平均 数 相 结合 的 位 置 更 新 策略 来 更 新 位 置 。 
1.4.1 BAR RAIR A 

a a dda 的 智能 算法 。 灰 狼 


种 群 有 严格 的 等 级 制度 ， 社 会 高 为 头 狼 a， 它 不 一 定 是 
最 强 的 狼 ， 但 必须 有 强大 的 py 负责 狩猎 、 捕 猎 、 食 物 


的 分 配 等 等 。 当 缺失 a IRIN, BA oa 狼 。 紧 接着 是 5 狼 ，o 
狼 追 随 ， 负 责 平衡 种 群 的 内 部 关系 ， 在 狩猎 时 实现 对 猎物 的 包 
AE 


灰 狼 狩猎 时 ， 根 据 收敛 因子 a 计算 系数 变量 A ， 当 |4|>1 
时 ， 灰 狼 的 群体 扩大 搜索 范围 ， 寻 找 更 好 的 食物 ， 这 对 应 于 全 
局 搜索 ;|4 |<1 时 ， 灰 狼 群 体 将 包围 圈 缩小 ， 这 对 应 于 局 部 搜 
Ro KAAF a 的 计算 如 下 所 示 ， 


“2 二 (13) 


其 中 :7 Ay OB, MAX 为 最 大 和 迭代 次 数 ，4 从 2 线性 递 
减 到 0。 搜索 猎物 时 ， 根 据 猫 物 位 置 X, 和 砍 狼 位 置 站 (7) 更 新 
灰 狼 个 体位 置 。 判 断 到 猎物 的 位 置 后 ， 根 据 a，B，5 狼 的 位 置 
Xo Xs, X Hi SHIM D, ，D; ，D; ， 对 猎物 
进行 包围 。 式 (14) ~ 16) 分 别 定义 了 o 狼 朝 a、B、5 狼 的 
前 进步 长 和 方向 。 


X,=X,-AD, (14) 
X,=X,—A,D, (15) 
X,=X,-A,D; (16) 


X (t+1) - Stoo Bs b 2x) 


(17) 
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a $, $: 


y 


灰 狼 搜索 过 程 中 ， 更 新 狼 个 体 的 位 置 


D=|CX,, (t)— X(t) (19) 

X (t+1)=X, (t)-AD (20) 

A=2ar -a (21) 

C=2r, (22) 

其 中 :C 和 4 为 系数 变量 ，D 为 灰 狼 个 体 与 食物 的 距离 ，X, 为 

MM, X(t) 为 第 t 代 时 灰 狼 个 体 的 位 置 ， 有 i 入 为 [0,1] 

的 随机 数 。 

2) 结合 调和 平均 数 和 静态 平均 数 的 位 置 更 新 

当 灰 狼 判断 到 猎物 的 大 概 位 置 时 ， 由 头 狼 a 带领 8B 狼 和 8 


狼 指导 狼 群 包围 猎物 ， 通 过 计算 a，B，5 狼 与 猎物 之 间 的 距离 
D., Ds, DRAWN, IID. 


)| (23) 


D, =|C,X, -xX (r) (24) 


D; =|C,X;-X (2) (25) 
其 中 : C,, C, CG AGA, X,, Xe, X 分 别 表 示 a. BA 
8 的 当前 位 置 。 根 据 式 〈14) - (16) 计算 @ RAO. BL 6 狼 的 
ESKAS X,, X,, X 
基本 GWO 算法 的 位 置 更 新 采用 静态 平均 数 方法 ， 即 取 适 
应 度 值 排名 前 三 的 个 体 进 行 位 置 的 平均 计算 ， 但 这 容易 导致 位 
置 更 新 不 灵活 ， 有 时 难以 找到 算法 的 最 优 解 。 为 此 本 文选 用 ] 了 
调和 平均 数 和 静态 平均 数 相 结合 的 方法 来 更 新 位 置 。 


3 
BO) aye x, (26) 


A (26) 和 静态 平均 数 ER AD ) 结合 完成 位 置 的 更 新 。 
设置 一 个 浆 值 ， 当 w 狼 和 8 狼 的 适应 度 差 值 大 于 立 值 时 ， 采 用 
调和 平均 数 更 新 位 置 ， 当 差 值 小 于 阔 值 时 ， 使 用 静态 平均 数 更 
新 位 置 。 经 反复 实验 验证 ， 本 文 阐 值 设置 为 0.1。 

1.4.3 IGWO 算法 实现 步骤 

a) 设 定 狼 群 变量 的 上 限 和 下 限 ， 在 上 下 限 中 随机 初始 化 狼 
群 位置 。 初 始 化 灰 狼 算法 的 参数 值 ， 包 括 种 群 规模 N, BKK 
代 次 数 max ，a、B、5 狼 的 位 置 ， 初 始 化 目标 函数 值 为 无 穷 

b) 计 算 每 个 狼 个 体 的 适应 度 值 ， 根 据 适应 度 值 的 大 小 来 确 


A (17) 定义 了 @ 狼 的 最 终 前 进 方向 。 当 4 递减 到 0 时 ， 灰 狼 

种 群 进行 攻击 ， 捕 获 猫 物 ，GWO 获得 最 优 解 。 

1.4.2 改进 灰 狼 算法 
1) 余弦 非 线 性 收敛 因子 


群体 智能 优化 算法 一 般 都 面临 着 全 局 搜索 和 局 部 搜索 之 间 
的 平衡 问题 ， 容 易 陷 入 局 部 最 优 。 由 式 (13) 可 知 ， 随 迭代 次 数 


的 增加 收敛 因子 4 由 2 线性 递减 到 0， 而 灰 狼 算法 在 收敛 过 程 
中 是 非 线 性 变化 的 ，a 的 线性 递减 不 能 完全 体现 算法 的 实际 搜 
索 过 程 。 于 是 ， 本 文采 用 了 余弦 非 线 性 收敛 因子 的 定义 方式 。 


azze af EY) (18) 


定 食物 的 位 置 。 选 择 适 应 度 值 排名 前 三 的 狼 个 体 作 为 g，B，5 
狠 。 

c) 搜 索 猫 物 时 ， 根 据 式 (20) 更 新 位 置 。 捕 获 过 程 中 ， 依 
据 式 (23) ~ (25) 更 新 a, B, ò 狼 的 位 置 。 根 据 式 (17) (26) 
Wa earner 8 狼 前 进 的 方向 和 距离 ， 更 新 整个 狼 群 的 位 


d) 保 存 迭 代 时 计算 的 最 优 个 体位 置 和 最 优 值 。 如 果 达 到 最 
大 和 返 代 次 数 ,输出 狼 群 的 最 优 个 体位 置 和 最 优 值 ;否则 , 返回 b)。 
1.5 改进 灰 a 

初始 化 灰 狼 种 群 ,每 个 灰 狼 个 体 的 位 置 由 要 优化 的 SVM 
参数 C 和 8 组 成 。C ee 即 对 误差 的 宽容 度 


录用 稿 


对 RBF 核 函 


数 的 参数 。 初 始 化 狼 群 的 数量 、 最 大 迭代 次 数 、C 


All g 取 值 的 范围 ,初始化 a、B、5 狼 的 位 置 和 狼 群 的 目标 函数 


值 。 
b) 将 初始 化 的 C 和 g 应 用 到 
将 测试 后 的 错误 率 作为 灰 狼 个 


为 


适应 度 值 / 


标 。 
c) 根 据 适 应 度 值 对 w、B、8 狼 的 目标 函 
函数 值 ， 则 将 a 狼 的 函数 值 更 新 为 适 


\F a 狠 的 目标 


到 SVM 的 样本 训练 和 测试 中 ， 
体 的 适应 度 值 ， 以 错误 率 最 小 化 


Wy 


应 度 值 ， 若 适 


则 将 5 狼 的 函数 值 更 
上 相应 灰 


狼 的 位 


d) 根 据 式 〈18) 更 新 收敛 因子 4 ， 
， 计 算 狼 群 位 置 更 新 的 步 长 和 方向 。 


中 。 


数值 进行 更 新 ， 若 


ee 数值 之 间 ， 则 将 B 狼 的 
目标 函 ens 应 度 值 ; 如 果 适 应 度 值 在 B 狼 和 8 狼 之 间 ， 


“a $, 
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8 GB, Windows10, MATLAB R2016b。 
2.1 PCMSFE 算法 的 性 能 分 析 实 验 


为 了 测试 PCMSFE 特 
个 采样 点 的 随机 噪声 和 1/f 噪声 ， 
特征 ， 其 中 的 参数 设 
相似 容 限 r =0.2*std (std 表示 时 
BE nr=2。 图 2 (a) ~ C) 分 别 表 
白 噪声 下 PCMSFE 和 


点 时 ， 随 机 


ERIH 


为 : 尺度 


能， 提取 512, 


Co) 分 别 表示 采样 点 为 512、 


PCMSFE 和 MSFE 对 比 结果 。 从 图 2 和 3 可 


为 适应 度 值 。 并 将 此 时 的 C 和 g 值 保存 


声 随 着 尺度 


复杂 度 趋 于 3 


F 稳 。 


的 增加 , 复杂 


HA 


BEAL 7 Al RAEL, 


计算 a、B、5 狼 的 位 
可 最 优 个 体 的 位 置 。 

e) 当 算法 达到 最 大 迭代 次 数 时 ， 输 出 a 狼 的 当前 位 置 和 最 
RE BI CE, g 值 和 训练 准确 率 。 


2 


2.2 


(Fi~Fo) 进行 仿真 实验 ，F1-Fo K 
Generalized Rosenbrock’s Function ， 


Generalized Rastrigin’s 


实验 结果 及 分 析 


多 尺度 模糊 业 


本 文 仿真 环境 为 Intel@ Core™ i5-8400M CPU @2.80 GHz, 


—$— PCMSFE 
一 是 一 MSFE 


N a 


L REO 
in 


= 


可 知 ， 在 随机 和 白 噪 


PCMSFE, 


PCMSFE, % 


MSFE 的 标准 
对 于 1/f 噪声 , KA 
线 更 


比 MSFE 


并 且 


随 着 


差 减 小 ,1 


和 1/f 噪声 
尺度 的 增加 ， 标 准 
个 尺度 上 产生 的 均值 误差 较 小 。 
点 增加 到 2048 点 时 ， 如 图 2 Cb) (ec) F 
总 体 上 MSFE 的 标准 


较 短 时 ， 采 


—$— PCMSFE 
—— MSFE 


0.5 


尺度 
(a) 采样 点 512 点 


图 2 随机 白 噪 j 


尺度 
(b) 采样 点 1024 点 


2 
还 
= a 
a 1.5 
BX 
K 
a 1 
—$— PCMSFE 六 —$— PCMSFE 
一 是 一 MSFE —— MSFE 
0.5 
10 15 5 10 15 


(a) 采样 点 512 点 


IGWO 算法 的 性 能 分 析 实验 


为 了 测试 IGWO 算法 的 性 能 , 选取 了 9 个 基准 


Function ， 


测试 函数 多 
数 分 别 为 Sphere Model， 
Quartic 


Function ， 


Function, Ackley’s Function, Six-Hump Camel-Back Function, 


Branin Function, Goldstein-Price Function。 其 中 Fi~F3 PK 


测试 函数 ，F4~Fs RK 


We 


EME 


Generalized Griewank 


Cb) 采样 点 1024 点 


BRILL 


AKA 
MSFE 对 比 结果 。 
1024, 2048 RIN, I/f 噪声 下 
以 看 出 ， 


ERAT 


1024, 2048 


来 对 比 MSFE 和 PCMSFE 
Alf t=1~15, WAH% m=2, 
间 J 列 的 标 ; A 


差 ) ， 边 界 梯 


点 为 512、1024、2048 


图 3 Ca) ~ 


ENLA NR 


度 逐 渐 降 低 , Lf 噪声 随 着 尺度 的 增加 ， 
采样 点 为 512 点 时 ， 由 图 2 (a) 和 图 3 Ca) 
H, MSFE 的 标准 差 大 于 
差 越 来 越 大 。 对 于 
当 采 样 点 


1024 


[图 3 (b) (c) 所 示 ， 


点 为 2048 点 时 , 如 图 3 Cc) 所 示 


加 平稳 ， 


示 准 差 差 别 不 大 。 由 此 说 明 数据 长 度 


| PCMSFE 算法 要 优 于 MSFE 算法 。 


—®— PCMSFE 
一 一- MSFE 


F PCMSFE 和 MSFE 对 比 


(c) 采样 点 2048 点 


图 3 1 噪声 下 PCMSFE 和 MSFE 对 比 


是 固定 维度 多 峰 基准 测试 


设置 为 ， 种 


本 文采 


IGWO1) , 


只 改进 位 


分 别 独立 进行 30 次 实验 , 记录 四 种 算法 的 


函数 是 单 
函数 是 多 峰 基 准 测试 函数 ，F7~Fo 函数 


éb 


价 算法 的 性 
峰 基 准 测 试 函数 中 ,IGWO1 在 


2 
€ 
z 
ak 
HX 
K 
ESN 1 
—$— PCMSFE 
一 里 一 MSFE 
0.5 
10 15 
尺度 


Co) 采样 点 2048 点 


函数 。 为 了 确保 结果 的 公 习 
单 规模 为 30， 最 大 迭代 次 数 为 500。 
用 基本 GWO， 只 改进 收敛 
更 新 的 GWO( 记 为 IGWO2) 和 


E 差 大 于 PCMSFE。 


~, PCMSFE 


Xx 


F 性 ， 所 有 


因子 的 GWO ( 记 为 


IGWO 


均值 和 标 ; 


在 F 中 标准 差 


EAL) 


Ż 


` ` 
平 
= 


能 ， 比 较 结果 如 表 1 所 示 。 从 表 1 可 以 看 出 ， 在 单 
Fl 中 均值 和 标准 
小 于 GWO; IGWO 在 Fi 中 均值 比 GWO 


F GWO, 
更 接近 


录用 稿 


于 0, 在 三 个 函数 中 标准 差 均 小 于 GWO。 对 于 多 峰 基 准 测试 函 
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数 ,IGWO 在 Fa, Fs 中 均值 和 标准 差 均 达到 了 最 优 解 0;IGWO1、 
IGWO2、IGWO 在 Fe 上 的 均值 均 小 于 GWO 的 均值 ， 且 标准 差 
均 为 0。 对 于 固定 维度 多 峰 基 准 测试 函数 ，GWO、IGWO1、 

IGWO2、IGWO 在 F?,Fe 上 均 达 到 了 最 优 解 , 且 IGWO1 和 IGWO 


在 Fy 标准 差 达到 0,IGWO2 和 IGWO 在 Fo 中 标准 差 小 于 GWO。 


此 可 见 ，IGWO 在 大 多 数 函 数 上 具有 一 定 的 优势 ， 既 可 以 提 


Fl 


高 函数 的 搜索 精度 ， 又 可 以 使 寻 优 结果 更 加 稳定 ， 综 合 性 能 较 


好 。 


为 了 能 够 进一步 说 明 IGWO 的 优越 性 ， 


本 文 给 出 了 Fi, Fa 


EA 


=| 


BUNS ues Hee, MA 4 所 示 。 从 图 4 可 以 看 出 ，IGWOI1 
以 显著 提高 GWO 的 收敛 速度 和 收敛 精度 ，IGWO 比 GWO 


上 共有 更 快 的 收敛 速度 和 更 高 的 收敛 精度 ， 说 明 本 文 提出 的 
IGWO 的 优化 算法 是 有 效 的 。 


1 


9 个 基准 测试 函数 的 测试 结果 


X 
GWO IGWO1 IGWO2 IGWO 
均值 标准 差 均值 标准 差 均值 标准 差 均值 标准 差 

Fi 1.80e-27 2.89e-27 4.49e-39 6.06e-39 1.43e-27 2.50e-27 1.58e-39 4.27e-39 
Fp 27.48 0.7615 26.7767 0.7593 27.0474 0.6543 26.9415 0.5963 
F; 0.0019 0.0012 0.0013 6.49e-04 0.0016 8.47e-04 0.0012 6.10e-04 
Fa 2.2455 9.03e-16 5.92e-17 3.24e-16 4.91e-11 0 0 0 

Fs 0.0246 0 0 0 0 0 0 0 

Fo 1.11e-13 0 1.46e-14 0 8.61e-14 0 1.27e-14 0 

F7 -1.0316 4.51e-16 -1.0316 0 -1.0316 4.51e-16 -1.0316 0 

Fs 0.3979 0 0.398 4.37e-04 0.3979 0 0.3979 0 

Fo 3 1.35e-15 3 1.35e-15 3 4.51e-16 3 4.51e-16 

010 10° 


目标 函数 值 


100 200 300 


迭代 次 数 


400 500 


2.3 ” 脑 电 信号 情感 识别 实验 


本 文 的 实验 样本 来 源 了 


图 4 ”函数 的 寻 优 收 敛 


F DEAP 数据 库 ,使 用 带 通 FIR 滤波 
人 信号， 去除 小 于 8Hz、 大 于 30Hz 的 噪声 。 
实验 分 为 效 价 、 唤 醒 度 、 优 势 度 、 喜 欢度 四 


—— GWO 


0 100 200 300 400 


TERME 


线 ( 左 为 Fl BAN F4 函数 ) 


我 们 将 


个 类 别 ， 同 时 ， 分 


别 比 较 在 高 / 低 喜 欢度 下 效 价 和 唤醒 度 的 情感 识别 率 。 因 DEAP 


数据 库 不 是 专门 用 来 评价 喜欢 度 的 数据 库 ， 有 些 被 试 者 的 数据 


中 喜欢 度 标 签 存 在 不 了 


每 个 被 试 者 


为 了 验 订 


统一 设 


i 
| 在 两 个 数 
取 总 数据 


F 衡 的 状况 , 于 是 本 文 忽略 了 高 / 
比率 小 于 0.3 的 被 试 者 的 数据 。 最 终 ， 有 21 个 被 试 者 的 数据 被 
用 于 测试 喜欢 度 对 效 价 和 唤醒 度 影响 的 实验 。 将 数 扩 
欢度 的 高 低 分 为 高 喜欢 度数 据 和 低 喜 欢度 数据 ， 分 另 
据 库 中 进行 效 价 和 唤醒 度 分 类 的 实验 ， 训 练 数 据 随 必 
的 80%, 剩余 的 20% 为 测试 数据 。 为 了 保证 实验 数据 的 公平 性 ， 
的 数据 以 标签 数 少 的 数据 量 为 标准 。 


L 


氏 喜 欢度 


库 按 照 喜 


EIGWO-SVM 模型 在 脑 电信 号 情感 识别 中 的 性 能 ， 
本 文 将 其 与 GWO、PSO、GA、FA 优化 的 SVM 模型 进行 比较 。 


优化 算 沪 


参数 设 


il 


500 


数 范围 设置 为 : CE[0.01,100],gE[0.01,100]。 通 过 寻 优 算法 训 
练 数据 集 ， 建 立 IGWO-SVM (图 表 中 标记 为 OURS) ， 
GWO-SVM, PSO-SVM, GA-SVM, FA-SVM, SVM 分 类 模型 。 


同时 , 本 文 建立 了 传统 的 MSFE 在 IGWO-SVM 分 类 下 的 模型 ， 


记 为 MSFE-IGWO 模型 。 


几 个 分 类 模型 的 平均 识别 率 对 比如 图 5 所 示 ，GWO-SVM 


模型 的 情感 识别 率 高 于 其 他 模型 的 情感 识 另 


对 SVM 优化 的 分 类 效果 高 


上 率 ， 说 明 灰 狼 算法 
F 其 他 各 个 算法 对 SVM 的 优化 效 


果 。IGWO-SVM 模型 在 效 价 、 唤 醒 度 、 优 势 度 和 喜欢 度 上 的 分 


类 效果 均 高 于 其 他 模型 , 证 实 本 文 对 于 GWO 的 改进 是 有 效 的 。 


MSFE-IGWO 模型 在 各 个 方 


的 种 群 规模 为 20， 最 大 迭代 次 数 为 200， 此 


E 了 寻 优 算法 都 可 以 达到 收敛 状态 SVM 的 寻 优 参 


IGWO-SVM 模型 ， 表 示 本 文 使 
征 更 具有 识别 性 。 表 2 为 每 个 模型 
K 2 可 知 ,IGWO-SVM 模型 识别 一 名 被 试 者 的 情感 的 时 间 


E 
FE 
最 短 的 ， 说 明 IGWO-SVM 模型 可 以 缩短 运行 时 间 ， 提 高 情感 


看 的 情感 识别 率 均 低 于 


的 PCMSFE 特征 比 MSFE 特 
对 一 位 被 试 者 识别 的 时 间 。 


m 


识别 的 效率 。 表 3 为 每 个 模型 在 效 价 、 唤 醒 度 、 优 势 度 、 喜 欢 


度 下 的 情感 识别 率 的 平均 值 、 标 准 差 、 最 大 值 和 最 小 值 。 从 表 


ChinaXiv 全 作 甘 明王 | 
. _ c NINAX IVA FA, JTY 
录用 稿 a 雪 ， 等 : DRLSS RAR MI IGWO-SVM 的 脑 电 情感 识别 


3 可 知 , IGWO-SVM 模型 在 各 个 分 类 上 的 平均 情感 识别 率 均 大 ” 低 对 效 价 和 唤醒 度 的 分 类 具有 一 定 的 影响 ， 观 看 低 喜欢 度 的 视 
于 其 他 模型 ， 标 准 差 值 小 于 其 他 大 多 数 模型 的 标准 差 ， 最 高 情 。 频 会 有 更 高 的 情感 识别 率 。 


感 识 别 率 可 达到 100%。 表 2 各 个 模型 识别 时 间 
喜欢 度 的 不 同 对 效 价 、 唤 醒 度 的 影响 如 图 6 所 示 ， 低 喜欢 模型 OURS GWO-SVMFA- SVM PSO- SVM GA- SVM 
度 下 的 效 价 和 唤醒 度 的 分 类 效果 比 高 喜欢 度 下 的 分 类 效果 好 ， 时 间 Cs) 6.778 7.403 7.123 14.646 35.862 
而 且 比 在 所 有 标签 下 的 情感 识别 率 高 。 可 以 看 出 ， 喜 欢度 的 高 
100 7 MB ce 
GE 高 喜欢 度 
四 所 有 标签 
90 F 
孙 
Og 80; 
jid 
{i 
70 
R 3 60 
效 价 ”唤醒 度 ”优势 ”喜欢 度 效 价 唤醒 度 
分 类 标准 分 类 标准 
图 5 四 个 分 类 标准 下 各 个 模型 情感 识别 率 对 比 图 6 高 / 低 喜欢 度 下 效 价 和 唤醒 度 分 类 
KI 各 个 模型 在 效 价 、 唤 醒 度 、 优 势 度 、 喜 欢度 下 的 分 类 效果 
分 类 标准 OURS GWO- SVM FA-SVM PSO- SVM GA- SVM SVM MSFE- IGWO 
平均 值 87.27 86.25 84.84 80.23 79.92 67.25 84.92 
seth 标准 差 5.832 5.923 6.155 6.879 7.580 5.880 5.626 
RUT 
最 大 值 100 97.5 100 92.5 95 81 92.5 
最 小 值 75 72.5 72.5 67.5 67.5 57 72.5 
平均 值 87.81 86.17 85.16 80.31 81.17 69.59 85.70 
标准 差 5.379 5.499 5.783 6.246 6.985 7.820 6.514 
唤醒 度 l 
最 大 值 100 100 97.5 90 95 85 97.5 
最 小 值 75 75 75 67.5 65 56 72.5 
平均 值 89.06 88.44 87.73 82.89 83.75 70.09 87.11 
标准 差 5.950 6.049 5.550 7.435 6.780 10.69 6.633 
优势 度 
最 大 值 100 100 100 100 100 100 100 
最 小 值 77.5 75 77.5 65 70 59.5 75 
平均 值 87.58 86.33 80.15 81.48 81.56 71.91 86.48 
标准 差 6.397 6.690 10.32 8.703 9.370 8.930 6.564 
喜欢 度 
最 大 值 97.5 97.5 95 97.5 95 95 97.5 
最 小 值 75 72.5 55 62.5 62.5 57.5 72.5 


试 的 情绪 状态 ， 提 出 的 算法 为 脑 电 信号 在 人 机 交互 、 辅 助 测 谎 
等 领域 的 应 用 打下 了 基础 。 
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